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Abstract 

The b — > 57 transition proceeds by a loop "penguin" diagram. It may be used to measure precisely 
the couplings of the top quark and to search for the effects of any new particles appearing in the 
loop. We present a preliminary measurement of the branching fraction of the exclusive decay, 
B o _> R *o 7 _ We uge g g x 1Q 6 B g decays to measure B{B° -» K*°j) = (5.4 ± 0.8 ± 0.5) x 10~ 5 . 
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Figure 1: The leading order Feynman diagram for the B° -» K*°-f decay 



In the Standard Model the exclusive decay B — ► i^* 7 proceeds by the b — > 57 loop "penguin" 
diagram shown in Fig. [I]. Precise measurements of decay modes involving these transitions and 
modes with the related b — ► c?7 transition such as -B — > /J7 will allow measurements of the 
top quark couplings Vt s and Vtd- The strength of these transitions may also be enhanced by 
the presence of non-Standard Model contributions JIJ. In minimal supersymmetric models, for 
example, the W can be replaced by a charged Higgs leading to significant enhancements in the 
branching fraction. In the first year of running the BABAR experiment has accumulated a dataset 
comparable to the world's largest to date, and this will increase by an order of magnitude over the 
next few years. The large dataset available and the state of the art detection systems employed 
at BABAR will allow measurements and searches of unprecedented precision. A comprehensive 
program to study these decays is now underway. The first step in this program is the preliminary 
measurement of the branching fraction of the exclusive decay mode B° — > K*°j using the leading 
decay mode, K*° — > K + ir~" . Here K*° refers to the if* (892) resonance, and charge conjugate 
channels are assumed throughout. The most precise measurement of the branching fraction to 
date, B(B° -» K*°~f) = (4.55lg;g ± 0.34) x 1(T 5 , is from the CLEO collaboration § and is in 
agreement with Standard Model predictions of (3.3 — 6.3) x 10~ 5 ||. 

The data were collected with the BABAR detector ||] at the PEP-II asymmetric e + e~ storage 
ring |5[. The results presented in this paper are based upon an integrated luminosity of 7.5 fb -1 of 
data corresponding to 8.6 x 10 6 BB meson pairs recorded at the T(4S) energy ( "on-resonance" ) 
and l.lfb -1 below the T{4S) energy ("off-resonance"). The BABAR detector simulation is based 
upon GEANT Q and tuned with data. Events taken from random triggers are used to measure the 
beam backgrounds. They are mixed into the simulated events, so that small changes in detector 
conditions, including dead channels, are also simulated. The simulated events are processed in the 
same manner as data. 

The selection criteria for this analysis have been optimized to maximize S 2 /(S + B) where S 
is the number of signal candidates expected, assuming B{B° -> K*°j) = 4.55 x 1(T 5 , and B is the 
expected number of background candidates determined from Monte Carlo. We compute quantities 
in both the laboratory frame and the rest frame of the T(4S). Quantities computed in the rest 
frame are denoted by an asterisk; e.g. is the energy of the e + and e~ beams which are symmetric 
in the T(4S) rest frame. 
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We begin the selection by requiring a high energy photon candidate in the calorimeter. The 
calorimeter consists of 6580 thallium doped Csl crystals arranged in a non-projective barrel and 
forward endcap geometry. We require the photon candidate to have an energy, measured in the 
laboratory frame, between 1.5 and 4.5 GeV and further require 2.20 < E* < 2.85 GeV. A photon 
candidate is defined as a calorimeter cluster [Q] in a region of good calorimetry (—0.74 < cos 9 < 
0.93), where 9 is the polar angle to the beam axis. The cluster must be isolated from any track. 
The distribution of individual crystal energies in the cluster is required to have only one maximum 
and a lateral profile consistent with a photon shower. These requirements remove backgrounds 
from high energy 7r° and r\ mesons where the two photons from the decay have merged into one 
cluster. In addition we form the combination of the high energy photon candidate with all other 
photons in the event with energy greater than 50 (250) MeV and require the invariant mass of the 
combination not to lie within two standard deviations of the known vr°(r/) mass. 

We next reconstruct the K*° from K + and tt~ candidates, by considering all pairs of tracks 
in the event. The tracks are required to be well reconstructed in the tracking detectors and to 
originate from a vertex consistent with the e + e~ interaction point. The tracking detector consists 
of a five layer double-sided silicon vertex tracker plus a forty layer stereo-axial drift chamber. A 
track is identified as a kaon if it is projected to pass through the fiducial volume of the particle 
identification detector. This detector is a novel ring imaging Cherenkov detector (DIRC) in which 
the radiating medium, an elongated quartz bar, also acts as a light guide to project the light 
cone onto an array of photo-multiplier tubes at the backward end of the detector. The high light 
yield and low mass of this detector are a significant advance in particle identification over previous 
experiments. We require a cone of Cherenkov light consistent in time and angle with a kaon of the 
measured track momentum. The K*° reconstruction is completed by requiring the invariant mass 
of the candidate pairs to be within 90 MeV/ c 2 of the K*° mass: 806 < M K+n - < 986 MeV/ c 2 . 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

cos(® T ) 



Figure 2: The event shape variable cos 9^ for B — » K* 7, K — > K + tt Monte Carlo and off 
resonance data. 

The B° candidates are reconstructed from the K*° and 7 candidates. There are backgrounds 
from continuum qq production with the high energy photon originating from initial state radiation 
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or from a 7r° or rj. We exploit event topology differences between signal and background to reduce 
the continuum contribution. The thrust vector T* of the event is computed, excluding the B° 
daughter candidates. Fig. Q shows the distribution of cos 9^ for signal Monte Carlo events and off 
resonance data, where 9^ is the angle between the high energy photon candidate and T*. In this 
frame the BB pairs are produced approximately at rest and therefore decay isotropically with a 
fiat distribution in cos 9^. The qq pair is produced above threshold recoiling against each other in a 
jet-like topology, which results in a cos6>^ distribution peaking at ±1 . We require | cos6^| < 0.78. 
Backgrounds are further suppressed by constructing two additional event shape variables which are 
uncorrelated to cos 9^. The angle of the B° candidate direction with respect to the beam axis, 9* B 
follows a 1 — cos 2 9b distribution and is flat for qq production. We require |cos6^| < 0.76. The 
helicity of the K*° decay, 9* H , is defined as the angle of the K + in the K*° rest frame with respect 
to the flight direction of the K*° in the T(4S) rest frame. This follows a 1 — cos 2 9* H distribution 
in B° — » K* 0/ y, K*° — ► K + ir~ , whereas the qq background is flat. We require | cos#J^| < 0.7. 

Since the B° mesons are produced via e + e~ — ► T(4S) — > BB, the energy of the B° is given by 
the beam energy, Ef. The beam energy is measured much more precisely than the energy of the B° 
candidate daughter particles. In particular for the decay B° — > K* ^y the energy resolution of the 
high energy photon dominates the measured B° candidate energy. We reconstruct the B° candidate 
substituting E% for the measured energy of the candidate daughters. We define the difference of 
the beam energy and energy of the B° daughters, AE* = E* K , + E* — E£. The B° mass is given 



by ities = Y^fe 2 ~~ \Pb \ i where \p* B \ is the momentum of the B° candidate calculated using the 
measured momenta of the charged daughters and the energy of the photon. In the calculation of 
mES we rescale the measured photon energy, E* by a factor k so that E£ = E* K , + kE*. This 
procedure corrects for the low energy tail in the E* distribution due to the incomplete containment 
of showers in the calorimeter. Using mES is an order of magnitude more precise than using the 
invariant mass. In addition the rescaling of E* enhances the resolution on mES by 20%. The 
resolution of 3.0MeV/c 2 is dominated by the beam energy spread. We select candidates with 
m E s > 5.2GeV/c 2 . 

Figure || shows the AE* versus w-es distribution for on-resonance data. Figure || shows the pro- 
jection onto the ttt-es axis, requiring —200 < AE* < 100 MeV. The ?tt,es distribution is fitted using 
an unbinned maximum likelihood fit. The background is empirically described by a function Q: 



dN / 1 wjes 2 



4-* 



where the parameter £ is measured from the off-resonance data sample which is required to pass 
the same selection criteria as the on-resonance data sample except that we remove the cos 9* H cut, 
and relax the selection requirements to | cos#J[ < 0.95 and |A£7*| < 500 MeV to gain statistics. We 
measure £ = 22 ± 13. The fit to the on-resonance data uses this fixed value of Q and adds a signal 
Gaussian whose mean and width are allowed to vary. We find a signal of 48.4 ± 7.3 events with the 
error coming from the statistical error of the fit. 

As a consistency check we plot in Fig. [| the AE* projection by requiring 5.274 < th-es < 
5.285 GeV/c 2 . The AE* distribution is fitted using a fixed shape determined from the Monte Carlo 
sample which is in good agreement with the data. We also plot M K + n - in Fig. ^ by requiring 
5.274 < m E s < 5.285 GeV/c 2 and -200 < AE* < 100 MeV. We fit using a Breit-Wigner shape and 
determine that the signal is consistent with coming from a K*°. 

The efficiency for the selection of B° — ► K* ^/ candidates is (15.5 ± 0.3)%. The branching 
fraction is determined using the yield, the efficiency and the total number of BB events in the 
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Figure 3: AE* versus m E s for B° -> K*°7,if*° -» K+vr- candidates from (8.6 ± 0.3) x 10 6 BB 
decays. The box indicates the boundaries used to make the projection of the AE*-m^ plane onto 
the mEs axis in Fig. |I[ 
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Figure 4: The m E s projection for B° -> K*°j, K*° -> K + ir- candidates from (8.6 ±0.3) x 10 6 BB 
decays. The fit uses a parameterized background function, described in the text, and a Gaussian 
for the signal. 



sample. The number of BB events is determined by counting the number of events passing a 
loose generic hadron requirement and subtracting the non-BB component estimated by scaling 
off-resonance data. We measure the number of BB events in the sample to be (8.6 db 0.3) x 10 6 . 
The branching fraction is measured to be B(B° — > K* ^) = (5.4 ± 0.8 ± 0.5) x 10~ 5 consistent 
both with previous measurements and with the Standard Model expectations. The first error is 
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Figure 5: The AE* projection for B° K*°j, K*° -> if +tt- candidates from (8.6 ± 0.3) x 10 6 
BB decays. 
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Figure 6: The M K+7T - distribution for B° K*°j , K*° -> K+vr- candidates from (8.6 ±0.3) x 10 6 
Si? decays. The fit is to a Breit-Wigner shape 



the statistical uncertainty in the yield. The second error is the systematic uncertainty. 

The total systematic error of 8.6% is a quadratic sum of several uncorrelated components 
tabulated in Table [lj. The systematic uncertainty in the background shape is obtained from varying 
the parameter £ within the allowed range from the off-resonance data. 

The primary check on track efficiency is obtained by studying the probability for observing drift 
chamber versus silicon detector-only tracks in inclusive D° — > K~ir + ir + ir~ . This is compared with 
the rate for finding the third track in one-versus-three topology tau-pair decays. A final check is 
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the observed multiplicity distribution in T(4S) events. A systematic error of 2.5% per track with 
Pt > 1 GeV/c on the overall efficiency scale is determined by comparing the three methods. In 
addition we measure the tracking resolution with e + e~ — ► p + p~ events and compare to Monte 
Carlo expectations. We find a difference in the tail of this distribution which leads to a systematic 
uncertainty in the K*° selection efficiency. 

The kaon identification efficiency in the DIRC is derived from a sample of D* + — > D°7r + in 
which the charge of the associated slow pion tags the flavor of the D meson and hence the charge 
of the kaon in the subsequent D — > K~ir + . The errors on this measured value are taken as a 
systematic uncertainty in the efficiency for kaon identification. 

The calorimeter energy resolution is measured in data using ir° and r\ meson decays, and e + e~ — > 
e + e _ 7 events. These are compared to Monte Carlo simulated events of the same processes. The 
differences result in a systematic uncertainty in the energy resolution which can broaden the AE* 
distribution causing an uncertainty in the efficiency of the AE* requirement. An overall energy scale 
uncertainty is estimated by using a data sample of i] meson decays with symmetric energy photons. 
The deviation in the reconstructed r] mass from the nominal i] mass estimates the uncertainty 
in the measured single photon energy. This resultant deviation in the central value of the AE* 
distribution causes a systematic uncertainty in the efficiency of the AE* requirement. The photon 
efficiency uncertainty is estimated using a sample of e + e~ — ► e + e~7 events compared to Monte 
Carlo simulated e + e~ — > e + e~7 events. The expected energy and position of the photon can be 
derived from a fit to the beam energy and the measured track momenta for the e + e~. We apply 
the photon selection cuts and compare the efficiency in the data and the Monte Carlo sample to 
derive the systematic uncertainty. The 7r°/r/ veto efficiency is tested by "embedding" a Monte 
Carlo generated photon into both off-resonance data and off-resonance Monte Carlo events. The 
difference in the measured efficiency is used as the systematic uncertainty. The merged 7r° modeling 
is tested by comparing a sample of merged 7r° in t~ — > p~v,p~ — > ir~ir° data and Monte Carlo. 
We assign a systematic in the photon selection efficiency by comparing the fraction of merged tt° 
removed in the two samples. 

Table 1: The fractional systematic uncertainties in the measurement of B(B° — ► K*°j). 



Uncertainty 


% of B(B" -> K^-i) 


Tracking efficiency 


5.0 


B counting 


3.6 


Kaon identification efficiency 


3.0 


Track resolution 


3.0 


Calorimeter energy resolution 


2.5 


Background shape 


2.3 


Monte Carlo Statistics 


1.9 


Calorimeter energy scale 


1.0 


Calorimeter efficiency 


1.0 


7T° /rj veto 


1.0 


Merged tt° modeling 


1.0 


Total 


8.6 



In conclusion, we report a preliminary measurement of the branching fraction of the rare decay 
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B(B° — > K*°7) = (5.4 ±0.8 ±0.5) x 10 5 with a precision comparable to previous experiments and 
consistent with Standard Model expectations. 
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